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摘要 : 碳 离子 束 治 疗 因 其 倒转 的 深度 剂量 分 布 和 较 高 的 相对 生物 学 效应 ,具有 


Fie", 


二 正常 组 织 损 


伤 小 、 副 作用 低 等 优势 ， 有望 成 为 头颈 部 肿瘤 放疗 的 主力 。 剂 量 算法 的 精确 性 直接 影响 临床 


照射 剂量 的 准确 性 ， 减 小 碳 离 子 束 治疗 剂量 计算 误差 有 助 于 提高 患者 生存 率 和 


局 部 控制 率 。 


本 文采 用 解析 笔 形 束 算法 对 5 例 头颈 部 患者 病例 进行 碳 离子 束 治疗 剂量 计算 ,3 


AUSF 


罗 模 拟 方法 对 治疗 计划 进行 剂量 验证 ,对 解析 计算 和 蒙特 卡 罗 模 拟 结果 进行 不 同 评价 标准 组 


合 下 (2%/2mm，3%/2mm，3%/3mm，3%/5mm) 的 伽 马 通过 率 分 析 、DVH 图 及 参考 指标 分 析 、 中 心 


轴 剂 量 分 布 对 比 。 结 果 显 示 : 笔 形 束 算法 在 3%/2mm 的 评价 标准 下 平均 伽 马 通过 率 


为 90. 441%, 


满足 碳 离 子 束 治疗 头颈 部 肿瘤 的 剂量 计算 要 求 。 笔 形 束 算法 在 束 流 横向 上 的 计算 结果 较为 精 


准 ,但 在 沿 束 流入 射 方向 上 的 计算 结果 受 路 径 上 组 织 非 均 匀 性 及 空 腔 的 影响 , 存 
为 降低 剂量 计算 误差 带 来 的 影响 , 碳 离子 束 治疗 头颈 部 肿瘤 时 需要 采用 更 保守 上 


在 射程 偏差。 
的 射 野 角度 选 


择 。 本 文 对 比 结果 可 用 于 指导 基于 解析 笔 形 束 剂量 算法 的 头颈 部 肿瘤 碳 离 子 束 治疗 计划 设 
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在 狭小 空间 内 集中 较 多 的 肌肉 、 骨 骼 、 血 管 和 神经 ， 各 种 重要 的 组 织 器 官 错综复杂 。 而 肿瘤 
大 多 形状 不 规则 , 与 周围 众多 重要 器 官 在 解剖 关系 上 不 可 避免 存在 重印 或 交叉 , 给 常规 手术 
治疗 带 来 严峻 挑战 中。 基于 生活 质量 和 功能 保护 的 考虑 ， 放 射 治 疗 是 头颈 部 肿瘤 治疗 的 主要 
手段 之 一 所。 近年 来 碳 离子 放射 治疗 (Carbon Ion Radiotherapy, CIRT) 因 其 良好 的 组 织 保护 和 
较 高 的 相对 生物 学 效应 在 头颈 部 肿瘤 患者 中 的 应 用 日 益 增 多 , 有 望 在 未 来 成 为 头颈 部 肿瘤 放 
疗 的 主力 B59。 

由 于 头颈 部 肿瘤 所 在 部 位 解剖 结构 特殊 ， 肿 瘤 周围 的 重要 器 官 受到 较 高 剂量 的 照射 会 
导致 严重 的 功能 损伤 。 因此 临床 治疗 计划 系统 (Treatment Planning System, TPS) 剂 量 计算 的 精 
确 性 极为 重要 ， 其 剂量 计算 结果 直接 影响 照射 剂量 的 准确 性 。 目 前 CIRT 主要 采用 解析 笔 形 
束 (Pencil Beam Algorithm, PBA) 算 法 进行 剂量 计算 , 其 计算 速度 较 快 , 计算 精度 也 相对 较 高 ， 
然而 其 不 能 很 好 地 处 理 非 均 匀 介 质 中 的 剂量 沉积 计算 外。 组 织 密度 高 梯度 变化 (如 空 腔 ) 的 
存在 会 导致 计算 结果 与 实验 测量 出 现 偏差 ,导致 浅 层 剂 量 以 及 射程 误差 四。 研究 表明 : 质子 
治疗 中 PBA 算法 在 头颈 、 肺 、 前 列 腺 等 具有 复杂 解剖 界面 部 位 的 剂量 计算 结果 有 一 定 误差 ， 
对 肿瘤 控制 率 的 影响 可 达 11%, 

以 往 对 于 PBA 算法 剂量 分 布 的 研究 对 象 主要 集中 为 水 模 体 、 非 均匀 平板 模 体 以 及 仿真 
模 体 ， 表 明 PBA 算法 在 非 均 匀 介 质 中 的 剂量 计算 结果 与 实际 测量 有 显著 性 差异 0。 考虑 
到 上 述 研究 中 仅 采 用 体 模 ， 难 以 反映 真实 患者 病例 的 复杂 情况 ,本 文 以 中 国 首 台 医用 重 离 
加 速 器 示范 装置 (型 号 : HIMM-01-GS-WW-01) 为 基础 ， 使 用 PBA 算法 对 5 例 头颈 部 肿瘤 患 
g 者 病例 进行 剂量 计算 ， 并 利用 被 用 作 剂量 验证 与 评估 的 金 标准 一 一 蒙特 卡 罗 (Monte Carlo, 
T MC) 模 拟 方法 对 治疗 计划 进行 了 剂量 验证 ,对 比分 析 PBA 算法 在 CIRT 头颈 部 肿瘤 的 剂量 计 

算 精确 性 ， 指 导 CIRT 临床 实践 中 的 治疗 计划 设计 。 
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2 材料 与 方法 

1 治疗 计划 软件 
本 文采 用 的 治疗 计划 软件 为 开源 MatRad 工具 包 ， 该 工具 包 使 用 MATLAB 脚本 编写 ， 
支持 光子 、 质 子 和 碳 离 子 的 调 强 放疗 治疗 计划 设计 。MatRad 支持 患者 DICOM 文件 的 导入 ， 
除了 勾画 器 官 外 , 可 以 实现 多 种 剂量 优化 算法 , 最 终 以 可 视 化 形式 生成 计划 评估 结果 与 报告 
033。 因 其 良好 的 可 读 性 与 易 改 性 ，MatRad 被 广泛 用 作 放 射 治疗 研究 的 测试 平台 。MatRad 
默认 使 用 的 解析 PBA 算法 是 质子 重 离子 领域 广泛 应 用 的 解析 算法 ， 其 基础 是 Fermi-Eyges 
理论 。 笔 形 束 算法 假设 照射 目标 是 由 不 同 介质 的 半 无 限 平板 层 层 堆 登 而 成 , 在 不 同 深度 处 笔 
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形 束 中 心 轴 周 围 为 均匀 介质 。MatRad 进行 粒子 


剂量 计算 时 需要 数据 库 中 不 同 能 量 水 中 剂量 


分 布 曲 线 数据 。 作 为 本 文 研 究 的 基础 ， 本 研究 团 


队 前 期 已 经 完成 基于 医用 重 离子 加 速 器 示范 


装置 (型 号 : HIMM-01-GS-WW-01) 点 扫描 束 流 线 对 应 的 数据 库 构 建 03]。 


2.2 蒙特 卡 罗 模 拟 工具 


本 文采 用 的 MC 软件 为 开源 GATE v9.1/GEANT4-11.0 工具 包 。 该 工具 包 基 于 


GEANT4(Geometry And Tracking) 内 核 开 发 , 其 主要 用 途 是 对 计算 机 断层 成 像 进行 MC 模拟 ， 


也 可 应 用 于 放射 治疗 领域 5491。 


(a) Nozzle 


(b) Ion Chamber (c) Ridge Filter 


ae HH- 


1 GATE 中 HIMM 装置 点 扫描 束 流 配送 系统 示意 图 


图 1 为 GATE 中 HIMM 点 扫描 束 流 配送 系 
流 强度 由 章 


统 模型 的 简要 示意 图 。 


流 穿 过 真空 窗 ， 束 


| 量 探 测 器 测量 , 束 流 的 形状 和 位 置 由 束 流 位 置 探 测 器 进行 测量 。 经 过 微型 硝 形 过 


滤器 后 ， 单 能 束 流 的 布拉格 峰 在 纵向 上 得 到 展 宽 ， 形 成 展 宽 为 3mm 的 微型 展 宽 布拉格 峰 。 


模拟 束 流 源 设置 如 表 1 所 示 。 


表 1 模拟 束 流 源 设置 
和 方向 立方 向 
能 量 扩散 束 斑 尺寸 义 方 向 束 流 发 射 度 立方 向 束 流 发 射 度 
束 流 角 散 度 束 流 角 散 度 
0.1% 3mm 3mrad 3mrad 10mm mrad 10mm mrad 


考虑 到 模拟 时 间 和 计算 精度 ， 仿 真 模拟 大 厅 、 柬 流 线 组 件 和 患者 模 体 之 间 的 伽 马 、 电 子 
和 正 电 子 产生 闵 值 分 别 设 置 为 In、lmm 和 0.1mm， 上 患者 模 体 中 的 步 长 限制 设置 为 0.1mm。 
本 文 使 用 的 物理 列表 为 : QGSP_BERT_ HP EMY, 其 中 包含 用 于 碳 离子 治疗 的 强 子 物理 过 程 


和 电磁 过 程 的 命令 。 模 拟 粒子 束 为 1 x 107, (EA 


服务 器 工作 站 上 的 CPU(Intel® Xeon® Gold 


6346 CPU@3.10 GHz, OS: Ubuntu 22.04 LTS) 进 行 计算 。 


2. 3 计划 设计 与 评估 方法 


随机 选取 甘肃 省 人 民 医 院 头 颈 部 肿瘤 患者 5 例 〈 数 据 来 源 ) ， 由 放疗 科 医 生 在 CT 图 
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部 计算 准确 性 。 因此 我 们 简化 了 计划 设计 的 流程 , 每 个 病例 仅 选 有 
表 1 为 患者 病例 的 具体 治疗 计划 参数 ， 所 有 的 治疗 计划 是 在 物理 剂量 的 基础 上 进行 优化 的 。 


究 的 主要 目的 是 对 比分 析 PBA 算法 在 解剖 界面 复杂 的 头颈 
区 1 NEK, 采用 单 野 照射 。 


束 斑 间距 和 治疗 床 角度 分 别 设置 为 3mm 和 0”， 设 置 30 个 治疗 分 次 。 
表 2 患者 病例 在 MatRad 中 具体 治疗 计划 参数 


病例 5 RE 感 兴趣 区 感 兴趣 区 种 类 剂量 目标 惩罚 因子 剂量 (Gy) 
PTV 

患者 1 270 PEX Squared Deviation 800 30 
CTV 
PGTV 

患者 2 60 PEX Squared Deviation 800 30 
GTV 
PGTV 

患者 3 90 PEX Squared Deviation 800 30 
GTV 
GTVp 

患者 4 270 P-GTV HE DX Squared Deviation 800 30 
P-GTV2500 
GTV 

患者 5 90 REK Squared Deviation 800 30 
PGTV 


(1) 1) Sy iE 28 2 NT 


每 个 患者 病例 在 MatRad 中 按照 相应 治 ; 


化 后 的 计划 数据 生成 对 应 的 GATE 粒子 源 


BH 
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将 MC 模拟 计算 结果 和 PBA 算法 计算 结果 归 一 化 处 理 


伽 马 通过 率 分 析 是 利用 数值 化 指标 训 


疗 计划 参数 进行 计划 设计 、 优 化 和 计算 。 根 据 优 
进行 MC 模拟 计算 069。 


后 进行 对 比分 析 研 究 。 


FE 估 两 种 剂量 分 布 之 间 的 一 致 性 ,通常 用 来 对 比分 析 
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对 局 部 或 全 局 最 大 剂量 点 归 一 化 的 剂量 差异 0729。 通 过 一 定 的 


当 y (7。) < 1， 则 该 剂量 点 通过 ，4da 


3%/5Smm) 进 行 全 局 剂量 的 伽 马 通过 率 


分 


测 剂量 点 与 基准 齐 


异 ) 限制 待 测 剂量 点 的 数量 ， 以 百分比 表示 的 伽 马 通过 率 代 表 两 种 剂量 
本 研究 采用 4 种 剂量 偏差 / 昌 


E. ` 
BAZ 


= 


TPS 剂量 计算 结果 和 实际 测量 结果 的 差异 071。y 通过 率 将 两 种 剂量 分 布 之 间 剂 量 差异 的 大 小 
Q 和 位 置 转化 为 统计 量 ， 其 定义 如 下 : 


Fw rp) = | 


(1) 


(2) 
间 的 距离 ，4D, 为 相 


评价 标准 


qu 


(剂量 差异 与 位 置 差 
分 布 之 间 的 一 致 性 。 


E 离 一 致 性 评价 标准 组 合 (2%/2mm, 3%/2mm, 3%/3mm, 


(2) 剂 量 体 积 直方 图 对 比 


选取 误区 作为 感 兴趣 区 域 ， 


Histogram, DVH) 并 比较 两 种 算法 


读 取 两 种 算法 计 香 


Dsow，Dos%，Dog%,) 的 剂量 差异 值 。 


(3) 剖 面 剂量 曲线 对 比 
提取 两 种 计算 结果 等 
果 差 异 。 


3 结果 与 讨论 
3.1 伽 马 通过 率 分 析 


对 两 种 计算 结果 分 别 使 月 


心 处 XYZ 轴 齐 面 剂量 分 布 ， 直 观 比较 PBA 算法 与 MC 模拟 结 


本 文 剂 量 差异 值 (%) 计 算 公 式 为 : pees 


马 通过 率 。5 例 患者 病例 使 


] PBA 算法 的 剂量 计算 结果 在 3%/5mm、3%/3mm 的 评价 


得 到 的 剂量 体积 直方 图 (Dose Volume 
F 均 剂量 (Dmwean) 和 不 同 百 分 比 目 标 体积 剂量 (Dw，Dsw， 


Pme X 100%. 


H 2%/2mm, 3%/2mm, 3%/3mm 和 3%/Smm 标准 来 评价 计划 个 


标准 下 


平均 伽 马 通过 率 均 之 95%, 在 3%/2mm 的 评价 标准 下 平均 伽 马 通 过 率 > 90%, 在 2%/2mm 的 


评价 标准 下 部 分 伽 马 通过 率 之 90%。 随 着 评价 标 


越 低 ， 详 见 表 3。 


表 3 5 例 患 者 病例 使 用 PBA FAM Mite 


准 越 来 越 严 格 , PBA 算法 的 验证 通过 率 越 来 


患者 编号 2%/2mm 3%/2mm 3%/3mm 3%/Smm 
1 90.236% 92.477% 94.490% 97.253% 

2 85.498% 89.993% 97.150% 99.085% 

3 84.591% 85.204% 92.103% 99.142% 

4 87.080% 88.782% 93.586% 98.401% 

3 94.088% 95.750% 98.920% 99.681% 
平均 值 88.299% 90.441% 95.250% 98.712% 


95% M>90% KA FE WAIT TX 


目前 对 常规 光子 调 强 放射 治疗 计划 ， 


[1920。 若 以 此 作为 参考 ，1 


到 碳 离 子 束 的 相对 生物 学 效应 ， 


的 可 
表 可 


90.441%， 最 优 效果 可 达 95.750%, PBA $ 


] 3%/2mm 的 评价 标准 ， 伽 马 通 过 率 应 > 
执行 性 ， 但 对 质子 重 离子 治疗 计划 未 提出 明确 的 评价 标准 
E 3%/2mm 的 评价 标准 下 平均 伽 马 通 过 率 为 


知 PBA 算法 


尤其 在 剂量 分 布 后 沿 布拉格 峰 末 端 具 
是 微小 的 距离 偏差 也 可 能 带 来 较 大 的 DVH 差别 。 因 此 对 于 碳 离子 束 的 
建议 采用 更 为 严格 的 评价 标准 。 


法 满足 CIRT 头颈 部 肿瘤 的 剂量 计 入 


要 求 。 考 虑 


高 的 LET (Ë 


， 即 使 


甸 床 治疗 剂量 验证 
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3.2 剂量 体积 直方 图 对 比分 析 
图 2 为 5 AAM RAA REKA DVH 图 , 其 中 实 线 代 表 MC 模拟 结果 ,虚线 代表 PBA 


算法 结果 。 表 4 AGRE DVH 统计 指标 的 相对 剂量 差异 值 。 除 4 号 患者 病例 外 ， 其 余 患 
者 病例 Dmean、D2、Ds、Dso 的 相对 剂量 差异 值 均 二 土 5%， 其 余 患 者 病例 Dos、Dog 的 相对 剂 
量 差异 值 大 部 分 过 土 5%;5 例 患者 病例 小 体积 受 量 (D。,，Ds)PBA 算法 计算 值 小 于 MC 模拟 值 ， 
大 体积 受 量 (Dos，Dog)PBA 算法 计算 值 大 于 MC 模拟 值 。 
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图 2 5 例 患 者 病例 靶 区 基于 MC 模拟 《〈 实 线 ) 与 PBA 算法 (虚线 ) 获得 的 DVH 图 


表 4 5 例 患 者 病例 PBA 算法 与 MC 模拟 结果 DVH 统计 指标 相对 误差 


患者 编号 K D mean Dz Ds Dso Dos Dog 
1 PTV 0.3% -3.3% -2.3% 0.8% 0.7% -2.0% 
CTV 1.0% -2.2% -2.0% 1.0% 4.1% 4.9% 
2 PGTV 2.0% -1.8% -1.0% 1.7% 7.2% 13.1% 
GTV 1.3% -1.3% -1.0% 1.4% 3.7% 4.2% 
3 GTV 3.0% -1.5% -0.8% 3.0% 71% 11.7% 
PGTV 3.3% -1.3% -0.7% 2.3% 7.1% 11.7% 
4 GTVp 7.3% -2.9% -1.9% 2.0% 83.2% 254.6% 
P-GTV 6.4% -3.4% -2.4% 1.4% 72.9% 249.6% 
P-GTV2500 7.6% -3.3% -2.2% 1.6% 105.5% 341.4% 
5 GTV -1.0% -5.1% -4.3% -0.5% 1.0% 1.2% 


PGTV -1.3% -5.0% -4.2% -1.0% 1.0% 5.9% 
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由 参考 指标 差异 值 可 知 : 对 大 部 分 病例 ，PBA 算法 和 MC 模拟 结果 对 同一 鞭 


考 指标 参数 不 同 ， 但 差异 范围 在 5% 以 内 ; PBA 算 


低估 ， 最 小 剂量 〈 即 大 体积 受 量 ) 存在 高 估 。 
量 误差 之 间 的 相关 诉 


为 基础 所 得 齐 


i 


E 不 明显 。 


为 全 局 伽 马 分 析 中 , 采用 全 局 最 高 剂量 进 


上 。 该 结果 也 说 明 单 独 使 用 伽 马 通过 率 进行 计划 验证 评估 剂 


3.3 XYZ 剖面 剂量 分 布 对 比 


结合 表 3 发 现 ， 全 


法 对 靶 区 最 大 剂量 〈 即 小 体积 受 量 ) 存在 


区 剂量 参 


局 伽 马 通过 率 与 DVH 参数 
如 4 号 患者 病例 在 2%/2mm 的 评价 标准 下 伽 马 
通过 率 为 87.080%， 明 显 优 于 部 分 患者 病例 ， 但 其 DVH 图 参考 指标 相对 差异 最 大 。 这 是 因 
行 归 一 化 , 其 分 析 结 果 不 能 反应 剂量 分 布 的 空间 信 
量 偏 差 是 不 严谨 的 。 


对 大 部 分 患者 病例 ， 两 种 结果 对 比 具 有 合理 的 剂量 一 致 性 。 图 3 为 剂量 一 致 性 较 差 的 4 
号 患者 病例 在 等 中 心 点 XYZ 轴 上 的 剂量 分 布 对 比 。 在 束 流 横向 方向 (YZ 轴 ) 上 ，PBA 算 


的 横向 半 影 剂量 计算 相当 准确 的 描述 相 


法 和 MC 模拟 剂量 分 布 具有 高 度 一 致 性 。 


BPN, ERIA H 


布 后 沿 位 置 ，PBA 算法 剂量 计算 结果 偏 


算 误 差 大 小 也 受 照射 路 径 复杂 度 的 影响 ,计算 误差 随 着 照射 路 径 上 器 官 组 织 、 


增加 。 


CIRT 的 剂量 学 优点 主要 是 具有 


端 布拉格 峰 区 。 


通过 治疗 计划 设计 ， 


照射 ， 而 正常 组 织 处 于 剂量 | 


Baines 


区 接受 低 剂量 照射 。 因 


X 


IE, CIRT 可 以 提供 精确 


[fini 


该 现象 也 与 文献 中 关于 PBA 算法 在 非 均匀 介质 中 
) ， 特 别 是 在 靶 区 剂量 分 
高 ， 明 显存 在 射程 偏差 。 如 图 4 所 示 ，PBA 剂量 
空 腔 的 增多 而 


ave 


到 转 的 剂量 分 布 且 射 程 固定 ， 大 部 分 剂量 沉积 在 射程 末 
肿瘤 靶 区 处 于 剂量 曲线 中 的 扩展 布拉格 峰 区 接受 高 剂 


pi 


MEK 


剂量 分 布 ， 同 时 保护 周围 正常 组 织 和 关键 器 官 。 然 而 PBA 算法 在 束 流 方向 的 剂量 计算 误差 


会 引入 碳 离子 束 射程 的 不 确定 性 ， 且 该 不 确定 性 


辑 区 剂量 覆盖 避免 潜在 风险 , 在 碳 离子 束 临床 治疗 中 


避免 照射 路 径 上 组 织 非 均匀 属 


变化 或 空 腔 ， 也 避免 关键 器 官 


Dose [Gy] 
o o o 
= 名 


0 150 2 
Y [mm] 


出 现在 束 流 剂量 分 布下 降 沿 。 


0 150 2 
Zimm] 


随 照 射 路 径 复 杂 度 的 提高 而 增 大 。 为 了 确保 
只 能 接受 更 保守 的 照射 角度 ,最 小 化 或 


图 3 4 号 患者 病例 (照射 角度 为 270”) 等 中 心 点 (a) X, (b) Y, (c) Z SHE MC 模拟 结果 ( 实 线 ) 与 


PBA 算法 (虚线 ) 计算 结果 剂量 分 布 对 比 


Dose 2 
Dose 2 


x [mm] 


100 200 
z[mm] 


z [mm] 


X 


4 不 同 患 者 照射 情况 对 比 。 (a) 1 号 患者 ， 低 复杂 度 照 射 路 径 ，(c) 4 号 患者 ， 高 复杂 度 照 射 路 径 


4 结论 

本 研究 基于 HIMM 装置 的 主动 点 扫描 束 流 配 送 方式 , 使 用 PBA 算法 进行 例 头颈 部 患 
者 的 碳 离子 放疗 治疗 计划 设计 ， 并 利用 MC 方法 对 治疗 计划 进行 了 剂量 验证 。 研 究 以 MC 
模拟 结果 为 基准 , 对 比分 析 两 种 结果 的 伽 马 通 过 率 、 DVH 图 及 参考 指标 和 XYZ 轴 剖面 剂量 
分 布 。 通 过 对 比分 析 ，PBA 算法 满足 头颈 部 肿瘤 的 CIRT 剂量 计算 基本 要 求 ， 其 在 束 流 横向 
上 计算 较为 精准 ， 但 在 束 流入 射 方向 上 的 剂量 分 布 存在 射程 偏差 ; 其 计算 结果 受 照射 路 径 上 
的 组 织 非 均匀 性 及 空 腔 影 响 。 由 于 头颈 部 肿瘤 的 结构 较为 复杂 , 为 降低 剂量 计算 误差 带 来 的 
影响 ， 进 行 CIRT 头颈 部 肿瘤 时 需要 采用 更 保守 的 照射 角度 ， 选 择 短 而 可 靠 的 照射 路 径 ， 最 
小 化 或 避免 照射 路 径 上 组 织 非 均匀 性 变化 或 空 腔 , 也 避免 危及 器 官 出 现在 束 流 剂量 分 布下 降 
沿 。 另 外 ， 通 过 金属 植 入 物 的 照射 路 径 也 应 该 尽 可 能 的 避免 。 


本 研究 通过 放射 治疗 剂量 计算 的 “ 金 标 准 ”MC 方法 对 实际 头颈 部 肿瘤 患者 病例 下 笔 形 
束 算法 剂量 计算 精确 性 进行 验证 ， 定 量 分 析 两 种 结果 的 差异 ， 可 用 于 指导 CIRT 头颈 部 肿瘤 
临床 上 的 治疗 计划 设计 。 
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Dosimetric Analysis of Analytic Pencil Beam Algorithm in 


Head and Neck Tumor Radiotherapy Plan Designing 


CHE Yuhang!”, HANG Shengcong!?, ZHANG Hui3“, LI Yazhou?**, LIU Xinguo*+°, FANG Tao!?, HU 
Qinyong!, LIQiang3 4 
( 1. Renmin Hospital of Wuhan Economic and Technological Development Zone (Hannan), Wuhan 430090, 
China; 
2. Wuhan University Heavy Ion Medicine Center, Wuhan 430090, China; 
3. Key Laboratory of Basic Research on Heavy Ion Radiation Application in Medicine, 
Gansu Province, Lanzhou 730000, China; 
4. University of Chinese Academy of Sciences, Beijing 100049, China; 


5. Gansu Provincial Hospital, Lanzhou 730000, China ) 
Abstract: Carbon ion radiotherapy is expected to become a workhorse for head and neck tumors radiotherapy due 
to its inverted depth dose distribution and high relative biological effectiveness, which has the advantages of less 
damage to normal tissues and low side effects. The accuracy of dose algorithms directly affects the accuracy of 
clinical irradiation doses, and reducing the dose calculation error of carbon ion radiotherapy could help improve 
patient’s survival and local control rates. In this paper, the dose distribution of carbon ion radiotherapy for 5 head 
and neck patients were calculated by the analytic pencil beam algorithm and validated using Monte Carlo method. 
Three analysis methods were carried out for the calculated and simulated results: (1) Gamma pass rate analysis 
under multiple evaluation criteria combination (2%/2mm, 3%/2mm, 3%/3mm, 3%/Smm), (2) DVH and reference 
index analysis, (3) central axis dose curve comparison. The analysis results showed that the average gamma pass 
rate of the pencil beam algorithm was 90.441% under the evaluation criteria of 3%/2mm, which met the dose 
calculation requirements of carbon ion radiotherapy for head and neck tumors. The pencil beam algorithm 
calculated more accurate in the beam transverse direction, and there is a range deviation by which the structure 
inhomogeneity and cavity exited in the irradiation path. In order to reduce the impact of dose calculation error, a 
more conservative angle selection to carbon-ion beam treatment of head and neck tumors is recommended. The 
results herein could be used to guide the design of carbon ion radiotherapy for head and neck tumors based on 


analytical pencil beam algorithm. 


Key words: Carbon ion radiotherapy, Pencil beam algorithm, Monte Carlo simulation, Head and neck tumor 
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